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The applicability, predictive power, and internal con- 
sistency of a modified BCS (MBCS) model suggested by 
Dang and Arima have been analyzed in details in [lj. 
It has been concluded that "The T-range of the MBCS 
applicability can be determined as far below the criti- 
cal temperature T c ", i.e., T « T c . Unfortunately, the 
source of our conclusions has been misrepresented in Q 
and referred to MBCS predictions at T » T c . 

Since above T c particles and holes contribute to an 
MBCS gap with opposite signs, the model results are 
rather sensitive to details of a single particle spectra 
(s.p.s.) (e.g., discussion in Sec. IV. A. 1. of 0). As so, 
it is indeed possible to find conditions when the MBCS 
simulates reasonable thermal behavior of a pairing gap. 
This can be achieved, e.g., by introducingsome particular 
T-dependence of the s.p.s. (entry 1 in [2j) or by adding 
an extra level to a picket fence model (PFM) (entry 2 in 
Q). But such results are very unstable and accordingly, 
the model has no predictive power. 

Dang and Arima explain poor MBCS results for the 
PFM (N = Q — 10) discussed in pj by referring to strong 
asymmetry in the line shape of the quasiparticle-numbcr 
fluctuations SAfj above T ~ 1.75 MeV (symmetry of SAfj 
is announced as a criterion of the MBCS applicability.) 
The space limitation is blamed for that in Q. Remember, 
particle-hole symmetry is an essential feature of the PFM 
with N = O. Thus, strong asymmetry is reported from 
the MBCS calculation in an ideally symmetric system. 

It has been found that a less symmetrical example 
N = 10, O = 11 satisfies better the MBCS criterion fy. 
Indeed, the model mimics behavior of a macroscopic the- 
ory in this case (see Fig. QJ)). But this example is the 
only one where the MBCS does not break, in a long row 
of physically very close examples with more limited or 
less limited s.p.s.. In all other examples we witness ei- 
ther negative heat capacity C v (Fig. QJi) or negative gap 
A (Fig. [IJ;) at rather moderate T (see also Q). 

Unfortunately, conclusion in 0] that "within extended 
configuration spaces . . . the MBCS is a good approxi- 
mation up to high T even for a system with N = 10 
particles" , is based on a single example while in all other 
N = 10 examples the MBCS yields unphysical predic- 
tions. 

The most serious problem of the MBCS is its thermo- 
dynamical inconsistency. It is not sufficient to declare 
two quantities < H >= Tr (HD) and £ representing the 
system energy, analytically equal by definition (as is done 
in footnote [8] of 0) to prove the model consistency. It 



> 



(a) 


(b) 


(c) 




— £1= 8 

— n= 9 
- ■-■ n = 10 


— £1=11 


— <\ 

"I 1 


- a = 12 

-0=13 

■ n = 14 - 






- \ 


































1 

1 
















'i i 







2 



1 2 1 

T (MeV) 



FIG. 1: (Color online) The MBCS pairing gap A (top panels) 
and specific heat C v (bottom panels) for the PFM with N = 
10 and (a) fi = 8,9, 10, (b) Q = 11, and (c) SI = 12, 13, 14. 
The pairing strength G = 0.4 MeV in all cases. 



is easy to find that expression for £mbcs (i n the form of 
Eq. (83) in (3) can be obtained in the same way as all 
other MBCS equations have been derived: straightfor- 
wardly replacing the Bogoliubov {uj, Vj} coefficients in 
£bcs(t=o) expression by {uj, Vj} coefficients. Numeric 
results in Fig. 9 of Q show that < H >mbcs and £mbcs 
have nothing in common while < H >bcs~ £bcs as it 
should be for thermodynamically consistent theory. 

Another example of the MBCS thermodynamical in- 
consistency is shown below. We calculate the system 
entropy S as 



Si 



dS 
~dt 



dt 



and 



S 2 = -£(2j + l) [i 



>,j hxrij + (1 — rij) In (1 



where rij are thermal quasiparticle occupation numbers. 
In Fig. |3 we compare Si and S2 quantities which refer 
to thermodynamical and statistical mechanical definition 
of entropy, respectively. The calculations have been per- 
formed for neutron system of 120 Sn with a realistic s.p.s. 

It is not possible to distinguish by eye Si and S2 in the 
FT-BCS calculation (solid curve in Fig.[21represents both 
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FIG. 2: Entropy of the neutron system in 120 Sn calculated 
within the FT-BCS (solid curve) and MBCS (dashed and dot- 
dashed curves). Notice the logarithmic y scale of the main 
figure and linear y scale of the insert. See text for details. 



quantities) as it should be for thermodynamically consis- 
tent theory. The MBCS Si and 5*2 quantities are shown 
by dashed and dot-dashed lines, respectively. They are 
different by orders of magnitude in the MBCS prediction. 

We stress that low T part is presented in Fig. [21 Dra- 
matic disagreement between Si (MBCS) and 5 2 (MBCS) 
representing the system entropy remains at higher T as 
well but we do not find it necessary to extend the plot: 
the model obviously does not describe correctly a heated 
system even at T ~ 200 keV. 

We show in the insert of Fig. another MBCS pre- 
diction: entropy Si decreases as temperature increases. 
This result is very stable against variation of the pair- 
ing strength G within a wide range and contradicts the 
second law of thermodynamics. 

Finally, we repeat, the conclusion in pj that "The T- 
range of the MBCS applicability can be determined as 
far below the critical temperature T c " is based on the 
analysis of the model predictions from T «T C and not 
on T >> T c results as presented in @]. 
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Forschungsgemeinschaft (SFB 634). 
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